XMM-Newton Observations Reveal Very High X-ray Luminosity 
from the Carbon-rich Wolf-Rayet Star WR 48a 

Svetozar A. Zhekov 1,4 , Marc Gagne 2 , and Stephen L. Skinner 3 



ABSTRACT 



We present XMM-Newton observations of the dusty Wolf-Rayet star WR 48a. 
This is the first detection of this object in X-rays. The XMM-Newton EPIC 
spectra are heavily absorbed and the presence of numerous strong emission lines 
indicates a thermal origin of the WR 48a X-ray emission, with dominant temper- 
ature components at kT coo ; pa 1 keV and kT/j 0i pa 3 keV , the hotter component 
dominating the observed flux. No significant X-ray variability was detected on 
time scales < 1 day. Although the distance to WR 48a is uncertain, if it is phys- 
ically associated with the open clusters Danks 1 and 2 at d ~ 4 kpc, then the 
resultant X-ray luminosity Lx ~ 10 35 ergs s _1 makes it the most X-ray luminous 
Wolf-Rayet star in the Galaxy detected so far, after the black-hole candidate Cyg 
X-3. We assume the following scenarios as the most likely explanation for the 
X-ray properties of WR 48a: (1) colliding stellar winds in a wide WR+O binary 
system, or in a hierarchical triple system with non-degenerate stellar components; 
(2) accretion shocks from the WR 48a wind onto a close companion (possibly a 
neutron star). More specific information about WR48a and its wind properties 
will be needed to distinguish between the above possibilities. 

Subject headings: stars: individual (WR 48a) — stars: Wolf-Rayet — X-rays: 
stars — shock waves 



Introduction 



As the early surveys revealed ( jPollockl 119871 ). Wolf-Rayet (WR) stars are sources of 
X-ray emission and massive WR+OB binaries are the brightest amongst them. Pointed 
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X-ray observations with modern observatories ( Chandra, XMM-Newton) conf i rmed this and 



provided us with grating; spectra of WR+OB binaries (ISkinner et al.l 1200 It iRaasen et al. 



20031 : ISchild et all 12004 ; IPollock et al.ll2005l : IZhekov fc Parkll2010bf ) that are rich in spectral 
lines from different ionic species which indicates a range of X-ray plasma temperatures. 
These findings suggest that the enha nced X-ray emission is pr oduc ed in colliding-stellar-w ind 
(CSW) shocks, as first proposed by iPrilutskii fc Usovl (119761 ) and ICherepashchuk I (119761 ). 



The available pointed Chandra and XMM-Newton observations suggest that the X-ray 
properties of presumably single WR stars might correlate with their subtype. Namely, the 
WN (nitrogen-rich) objects have X-ray spectra with prominent emission lin es originating 



from an admixture of cool (kT < 1 keV) and hot (kT > 2 keV) plasma (ISkinner et al. 
20101 ). Only one of the four WO (oxygen-rich) stars in the Galax y, WR 142, has been 



observed so far a nd it is a weak but extremely hard X-ray source (lOskinova et al.l 12009 



Sokal et al.ll2010l ). In contrast, the single WC (carbon-rich) stars are conspicuously X-ray 
quiet: all the Cha ndra and XMM-Newton pointed observat i ons of WCs so far h ave resulted 
in non-detections (lOskinova et al.ll2003r ISkinner et al .1120061 ; ISkinner et al.ll2009l ). The X-ray 
production mechanism in single WR stars is still not well-understood, but the ubiquitous 
presence of a high-temperature X-ray component in WR spectra is a clear i ndication that 



other processes besides embedded, radiative wind shocks are contributing ( ISkinner et al. 
201CH and references therein). 



Overall, binary systems, especially wide binaries, are more X-ray luminous than com- 



para ble single WR stars . For single WR stars Lx is typically in the range 10 31 — 10 33 ergs 



s _1 (ISkinner et al.ll2010l ). while for the known wide WR+OB binaries it is about an order of 
mag nitude higher, with WR 140 having the largest known luminosity of (1 — 2) x 10 34 ergs 
s _1 ( Pollock et al.ll2005l ). As a rule, if a WR star is detected with L x > 10 33 ergs s _1 , it is 
most likely a wide CSW binary system. 

In this Letter, we report the first X-ray detection of the Wolf-Rayet star WR 48a , which 
establishes it as a very luminous X-ray source and as such, a likely multiple system. 



2. The Wolf-Rayet Star WR 48a 



WR 48a was discov ered in a near-infrared survey and originally cl assified as a WC9 ob - 
ject (IDanks et al.lll983l ) with a current spectral classification of WC8 (Ivan der Huchtll200ll ). 
WR 48a is located inside the G305 star-forming region in the Scutum Crux arm of the Galaxy. 
Within 2' of WR 48a are two compact infrared clusters (Danks 1 and 2) and its proximity t o 
them suggests that this WR star likely originates from one or the other (IDanks et al.lll984j ). 
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The optical extinction towards WR 48a is very high, Ay = 9.2 mag (IDanks et al.lll983f). w ith 
a small fraction of it (~ 2 mag) coming from circumstellar material (IBaume et al.l 120091 ) . 



Th e infrared vari ability of WR 48a suggests that it is a long-period episodic dust- 
maker ([Williams! Il995l ) . but the dust-formation h i story indicates that the actual situation 
might be more complex. Namely, IWilliams et al.l (120031 ) proposed the existence of a short- 
period (P « 1 yr) binary in addition to the longer-period binary (P > 23 yr), where the 
latter is thought to cause the episodic dust formation. High-angul ar resolution near-in frared 
observations revealed a 70 mas size for the dust shell in WR 48a (jMonier et al.l 120071 ). 



The distance to WR 48a is not well constrained . If it is physically associated with 



the open clust e rs Da nks 1 and 2, IDanks et al.l (119831 ) estimated a distance of 4 kpc, but 
van der Huchtl (120011 ) adopted a photometric distance of 1.21 kpc based on a WC8 spectral 
classification. 

No information is available for the mass-loss rate and velocity of the stellar wind in 
WR 48a. But, indirect evidence of its powerful wind com es from mid-infrared observations 
of the region surrouding WR 48a (jClark fc Porter 1120041 ). which reveal a three- lobed neb- 
ula whose morphology is suggestive of a wind-blown structure and triggered massive star 
formation. 



3. Observations and Data Reduction 

WR 48a was observed with XMM-Newton on Jan 9, 2008. As discussed below, the 
presence of numerous spectral lines and our spectral analysis show that the X-ray emission 
is thermal and is heavily absorbed. The high absorption masks most emission in the RGS 
energy range below 2 keV, so we were not able to obtain any valuable information from the 
grating spectra. Thus, this study is based on the XMM-Newton European Photon Imaging 
Camera (EPIC) data which after excluding background flares resulted in effective exposures 
of 57 ksec for the pn spectrum and 71 ksec per MOS spectrum. The corresponding total 
number of X-ray counts in the (0.5 - 10 keV) energy range was: 75,000 (pn), 51,000 (MOS1) 
and 50,000 (MOS2). We used the XMM-Newton Science Analysis System (version 10.0.0) 
to extract the source and background spectra and to construct the corresponding response 
matrices and ancil lary response fi les. We used standard as well as custom models in version 



11.3.2 of XSPEC (lArnaud Ill996l ) for our spectral analysis. 



Figure [T] shows the soft and hard-band images of WR 48a that illustrate an excellent 
correspondence between the X-ray source and the optical coordinates of this star (J2000): 
a X MM = 13H2 m 39!47, 8 XMM = -62°42'56"00; a SIMBAD = 13 h 12 m 39!65, 5 SIMBAD = 
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— 62°42 / 55 / . / 80. There is an optical star located ~ 9" from WR 48a in the southwest di- 
rection (IWallace et al.l 120031 ) but we find no detectable X-ray emission from that source. 
Thus, all the EPIC X-ray emission should be associated with WR 48a0. 



4. Global Spectral Fits 



A general property of the XMM-Newton spectrum of WR 48a is that it is heavily 
absorbed and most of the X-ray photons are at energies > 1 keV. We thus re-binned the 
pn, MOS1, and MOS2 EPIC spectra in the 0.5 - 10 keV range to a minimum of 100 counts 
per bin and fitted them simultaneously. We used discrete-temperature plasma models to 
derive the global properties of the X-ray emitting plasma in WR 48a. We considered two 
limiting cases: (i) models that assume plasma in collisional ionization equilibrium (CIE), 
and (ii) plasma with non-equilibrium ionization (NEI). For consistency in the atomic data, 
we made use of emission models vapec and vpshock in XSPEC. Since no specific information 
is available on the WR 48a abundances, we adopted the following set of abundances typical 
for the WC stars (by number): H = 0.00, He = 0.618, C = 0.248, N = 0.00, O = 0.120, 
N e = 1.15 x 10~ 2 , Mg = 1 .68 x 10~ 3 , Si = 4.23 x 10~ 4 , S = 9.40 x 10" 5 , Fe = 2.36 x 10~ 4 



( Ivan der Hucht et al.l 119861 ) . Ar and Ca are not present in the Ivan der Hucht et al.l (119861 ) 
abundance set, so assumed a fiducial value of 1.2 x 10~ 5 . We varied the Ne, Mg, Si, S, Ar, 
Ca and Fe abundances in the spectral fits since these elements provide most of the X-rays at 
energies > 1 keV. Tabled] and Figure [2] present the results from the two-temperature model 
spectral fits of WR 48a. 

Though the quality of the two-temperature vapec model fit is accep t able, the NEI vp- 
shock model provides a better fit to the data. Applying the iGorensteinl ( 119751 ) conversion 



(N/f = 2.22 x 10 21 Ay cm -2 ), we find that the X-ray absorption is consistent (within 10-15%) 
with the optical extinction (Ay = 9.2 mag; § [2]). But, the X-ray data may indicate some 
extra absorption (~ 40 — 50%) if a more recent con yersion is used, ba s ed on studies of star - 



forming regions: N# = (1.6 — 1.7) x 10 Ay cm 2 (IVuong et al.ll2003l : iGetman et al.l 120051 ). 



The derived abundance values are in general consistent with the adopted set of non-solar WC 
abundances. However, the abundance values derived from X-ray analysis and especially from 
CCD spectra should be considered with caution. Also, due to the high X-ray absorption no 
valuable information could be derived for the abundances of carbon, nitrogen and oxygen. 



1 WR 48a is also detected with Chandra (Obsld: 8922) but pileup in the ACIS CCD is so high that no 
reliable spectral or timing information for WR 48a could be extracted. The ACIS image confirms that the 
WR 48a identification is correct (J2000: a chandra = 13' l 12 m 39!63, S C handra = -62°42'55'/90) and no other 
sources contribute to its X-ray emission (Gagne et al. 2010, in preparation). 
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But perhaps the most important result from the global fits is that at least two components 
with different plasma temperatures are needed to obtain an acceptable representation of the 
observed spectrum, and the hot component requires a relatively high temperature. The hot 
component (kT « 2.8 — 3.2 keV) provides about 75-85% of the observed flux, while the cooler 
component (kT 0.9 — 1.1 keV) provides ~ 60% of the unabsorbed X-ray emission from 
this WC star. 

In addition to the two-component models, we ran some more complex models: (i) a 
model with a thermal CIE plasma distribution based on our custom model which is similar 
to c6pvmkl in XSPEC, but uses the apec collisional plasma for the X-ray spectrum at a 
given temperature; and (ii) a model with a distribution of adiabatic NEI shocks which made 
use of our custom XSPEC m odel which was suc cessfully used in the analysis of the X-ray 



spectra of SNR 1987A (e.g., IZhekov et al.l 120091 and the references therein). Due to the 



high X-ray absorption, the X-ray spectral fits of WR 48a are not sensitive to the presence 
of very cool plasma. Thus, we considered only plasma temperatures above 0.5 keV for the 
emission measure distribution in WR 48a. The fit results are given in Table [TJ Figure [3] 
shows that the derived distribution of emission measure is bimodal, with temperature peaks 
closely matching the components of the two-temperature models. 

From the above results, it seems reasonable to conclude that there is evidence for a 
two-temperature distribution in the X-ray emitting region of WR 48a. But, it is worthwhile 
to treat such results with caution when based on CCD spectra since the spectral lines are not 
well resolved and the lines (and their ratios) are most sensitive to the presence of a variety of 
plasma temperatures. But, we believe that the bimodal distribution is in fact an indicator 
of a temperature-stratified X-ray emission region in WR 48a. 



Discussion 



The two most important results from the model fits to the XMM-Newton spectra of 
WR 48a are that it is a very luminous X-ray source and thermal plasmas with high temper- 
ature dominate its emission (Tableland Fig. [3]). Using the unabsorbed flux values from 
Table EJ the X-ray luminosity of WR 48a is L x (0.5 - 10 keV) = (0.5 - 2.1) x I0 34 d 2 kpc ergs 
s -1 , where dk pc is the distance in units of kpc. The upper limit is from the NEI shock model 
which gives a better quality of fit. The distance to WR 48a is not tightly constrained and a 
range of values is available in the literature: dk pc = 1.21—4 (§H]). But, based on its proximity 



to the open clusters Danks 1 and 2 and the sim ilar interstellar extinction (e.g.. iDanks et al. 



19831 : IClark fc Porter 1120041 ; iBaume et al.ll2009l ). the distance of 4 kpc seems more realistic 



which puts the X-ray luminosity in the range Lx = (0.8 — 3.4) x 10 ergs s . In this case 
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and excluding the most X-ray luminous WR star, Cyg X-3, which is a WR binary with a 
compact companion (a neutron star or a black hole), WR 48a is the most luminous WR star 
in the Galaxy among those observed so far. For exam ple, the X-ray luminosity of the bright- 



est C SW binary, WR 140, is: (0.5 - 3) x 10 34 ergs s~ l (IZhekov fc S 



kmner 



2000; Pollock et al. 



20051 ) . Adopting the bolometric luminosity from IClark fc Porter I (120041 ) , we note the very 
high value of IgLx/Lboi = [—4.3, —3.7] for WR 48a. All this raises the interesting question 
about the X-ray production mechanism in WR 48a. 

Colliding S tellar Winds. As already mentioned, WR 48a is an episodic dust-maker 
(jWiUiamslll995l ) which suggests that CSWs in a wide binary system with an eccentric orbit 
might provide most of its X-ra y emission, as is the c ase for the prototype episodic dust-maker 
amongst WR stars, WR 140 (jWilliams et al.lll990l ). Unfortunately, the lack of information 
about the stellar wind parameters of WR 48a did not allow us to make a comparison between 
the theoretical predictions based on hydrodynamic CSW modeling and observations. On the 
other hand, some simple (e.g. qualitative) considerations are possible. The temperature of 
the hot plasma component deduced from the global sp ectral fits (SHJ) can provide an estimate 
of the stellar wind velocity (see §5.2 in IZhekovl 120071 for discussion of CSW models versus 
discrete-temperature models ). For the case of typical WC abundances, the postshock plasma 
temperature is kT = 3.09 Vi 000 keV, where Viooo is the shock velocity in units of 1000 km 
s _1 . From Tableland Fig. [3j we see that the stellar wind of the WC star in WR 48a must 
have a velocity V W i n d > 1000 km s _1 . Such high velocities w ould be expected f or WC8-9 
stars since they have typical a verage wind speeds 1400 km s _1 (IPrinja et al.lll990f ). 1300 km 



s 1 ( lEenens fc Williams! Il994[ ). We note that the derived values for ionization age of the 



shocks in the two-temperature model (Table [T]) and in the model with a distribution of NEI 
shocks (Fig. [3]) are also qualitatively consistent with the CSW picture. Namely, the higher 
density and higher temperature plasma is located near the axis of symmetry (the line-of- 
centers between the two stars) while the less dense and cooler plasma is found downstream 
from that axis. Thus, it is natural in the CSW picture that the higher temperature plasma 
will have a larger ionization age. We can use the scaling law for the CSW X-ray luminosity 
with the mass- loss rate (M), wind velocity ( v) and binary separation (D): Lx oc M 2 v~ 3 D~ 1 
( ILuo et al.lll990l ; iMyasnikov fc Zhekovlll993l ) for a comparison between WR 48a and the clas- 



sical CSW binary W R 140. Adopting the WR 140 binary parameters ( [Williams et al.lll990 



Pollock et all 120051 ). a wind velocity of 1300 - 1400 km s" 1 in WR 48a and assuming that 
its mass loss is equal to that of WR 140, we see that the binary separation in WR 48a could 
be similar to that in WR 140 even when the WR 48a X-ray luminosity is about an order of 
magnitude larger than that of WR 140 (see above). Therefore, the CSW picture in a wide 
binary system is at least qualitatively consistent with the observational data for WR 48a. 
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Magnetically Confined Wind Shocks (MCWS). This mechanism is capable of producing 
hard X-ray emission. It was proposed to explain the high plasma temperature found in X-ray 
emiss ion from young massive stars and req uires the presence of a relatively strong magnetic 
field (jSchulz et al.ll2003llGagne et al.ll2005l and references therein). Then, could it be that the 
WC star in WR 48a is a magnetized object or has a close magnetized companion? We note 
that the X-ray luminosity of the p rototype MCWS object, 9 1 Ori C, is Lx (0.5 — 10 keV) ~ 



10 33 ergs s 1 (IGagne et al.l 120051 ) which is one to two orders of magnitude below that of 



WR 48a. Thus, we can likely rule out the case of a close magnetized companion. On the 
other hand, the MCWS mechanism does not seem promising if adopted directly to the WC 
star in WR 48a. First, in order to be efficient this mechanism requires a relatively strong 
global magnetic field to confine the massive WR wind and no such fields have been reported at 
present for these stars. Second, because of the expected decay of the magnetic field strength 
with the age of a mas sive star, the MCW S mechanism is associated only with young massive 
stars (age < 1 Myr; ISchulz et al.l 120031 ). Thus, the likely association of WR 48a with the 
open clusters Danks 1, 2 and the recent estimate of ~ 5 Myr for their age (IBaume et al. 
20091 ) present serious difficulties for the MCWS model in the case of WR 48a. 



Wind Accretion Shocks. The presence of a close degenerate companion (e.g. a neutron 
star) seems very intrigui ng and accretion onto such an object can in general provide high 



X-ray luminosities (e.g., iDavidson &: Ostrikerl Il973f ). Details depend on the actual binary 



parameters and wind properties of the main stellar component which are yet unknown for 
WR 48a. But, this case may not be unlikely if WR 48a is indeed a member of the 5-Myr- 
old open clusters Danks 1, 2 and if it once had a more massive close companion which has 
already evolved and exploded as a supernova. 

We have listed above some physical mechanisms that might be able to explain the 
observed X-ray properties of WR 48a. One could continue this list by considering other 
mechanisms or even more complex combinations of those mentioned above. But, all such 
mechanisms and the corresponding physical picture remain quite speculative because of 
the scarcity of detailed observational data for WR 48a . Thus, we will end our discussion 
by simply mentioning the cases that look least speculative to us and briefly discuss what 
observational information may help distinguish between them. 

For the moment, the following cases are the most likely explanation for the X-ray prop- 
erties of WR 48a: (i) CSWs in a wide binary system; or (ii) accretion wind shocks in a 
close binary: WR+compact companion (a neutron star). Variability is a key parameter for 
distinguishing between these cases. Namely, if long-term (months, years) X-ray variability 
is established from future observations, then case (i) would be strongly favored. Further- 
more, case (i) would be strengthened if variable non-thermal radio (NTR) emission from 
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WR 48a is eventually discovered and even more so if the X-ray emission is found to modu- 
late at a similar period. In general, N TR WRs are commonly associated with wide binary 
systems ([Dougherty fc Williams! 120001 ). Alternatively, if short-term (a few days) variability 
is established or a sudden change of the X-ray luminosity is detected, case (ii) would be 
favored. The XMM-Newton data show no significant variability on timescales < 1 day; the 
Kolmogorov-Smirnov test rules out variability at t he 95% confidence lev el. On the other 
hand, based on the infrared variability of WR 48a, IWilliams et al.l ( 120031 ) proposed that it 
is a triple stellar system. If so, we may expect some very long-period variability due to 
the X-rays from CSWs in a wide binary with eccentric orbit and some with a much shorter 
period that could result from the stellar wind shocking onto a non-degenerate companion. 
In this case, WR 48a would be very similar to the CSW binary WR 147 which was recently 
resolved into a double X-ray so urce with Chandra, of w hich one component is variable and 
possibly an unresolved binary (IZhekov fc Park! l2010al .b). However, WR 48a has a much 
higher X-ray luminosity. Finally, we hope that the unusually high Lx of WR 48a will moti- 
vate deeper follow-up observations across the entire spectrum (radio, infrared, optical, UV, 
X-rays). Such observations would provide us with valuable information about the star and 
its wind properties that could be very helpful for constraining the physical picture in this 
remarkable but understudied Wolf-Rayet system. 
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Fig. 1. — The EPIC-pn images of WR 48a (linear intensity scale). R.A. (J2000) and decl. 
(J2000) are on the horizontal and vertical axes, respectively. The circled plus sign gives the 
optical position of WR 48a (SIMBAD). 
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Fig. 2. — XMM-Newton background-subtracted spectra of WR 48a and the two-temperature 
fit with NEI shocks (Tabled]). The MOS1 and MOS2 spectra are drawn in red and green, 
respectively. The low and high-temperature components are shown in the upper panel with 
dashed and dotted lines, respectively. 
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Fig. 3. — Emission measure (EM = J n e nH e dV) of the WR 48a distribution of thermal 
CEI plasma and of adiabatic NEI shocks for a reference distance of d = 1 kpc. The arrows 
indicate the plasma temperature from the two-temperature fits (see Tabled]). The dotted 
line shows the derived ionization age in the case of NEI shocks (n e t). 
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Table 1. Global Spectral Model Results 



Parameter 


2T vapec 


2T vpshock 


CIE Plasma 


NEI Shocks 


y 2 /dof 


1 382/1066 


1169/1064 


1403/1064 


1 173/1062 


N„ fin 22 rm - 2 ) 


2 27+2'2 3 , 

^•^ ' -0.02 


2 30 + 2'25 


^•^-0.03 


2 29+2 

z " z,a -0.04 


kTi (keV) 


87+2^ 

u -° ' -0.01 


i - uu -0.07 






kT 9 fkeV) 


3 26+S-2S 


9 09+O.O7 
z, - OZl — 0.06 






EMi (10 54 cm" 3 ) 


3 58+2^ 


2 42+°-?« 

^•^-0.16 






EMo (10 54 cm -3 ) 


4 21+S-S5 

4: - z '- L -0.04 








Ti (10 11 cm" 3 s) 




9 ^9+0.56 
z -^ z -0.42 






r 2 (10 11 cm" 3 s) 




8.09 + {; 3 ^ 






Ne 


u -^ u -0.09 






11+ - 03 

u - ±i -0.03 


Mg 


n no+0.04 
u - zo -0.03 


1 3 +0 - 02 

u - io -0.02 


25+ - 04 

u - zd -0.04 


1 3+ 002 

u - io -0.02 


Si 


i oc+0.06 
± - od -0.06 


64 +0 - 06 

u - u ^-0.05 


1 OO+0.07 
± -°°-0.07 


65+ - 06 

U - DO -0.06 


s 




1 78+0.15 


A OO+0.20 

4 -^-0.20 


1 01+0.12 

1.»1_q.12 


Ar 


Q 79+O.94 
a - ' z -0.92 


9 70+O.I9 
z - 'O-0.36 


y - do -0.94 


9 09+0. 38 
Z - OZ -0.34 


Ca 


K i q+0.65 
°- L J -0.62 


1 O7+0.34 
■ L - 3 ' -0.33 


r 99+O.69 
d - zz -0.70 


9 no +0 - 33 

z - uu -0.41 


Fe 


i o fi +0.06 


1 01 +0.07 
i - oi -0.07 


1 07+O.O7 
J-- ' -0.07 


1 ^o +0 - 07 

i - ou -0.07 


Fx (10 -11 ergs cm -2 s^ 1 ) 


1.06 (4.21) 


1.05 (17.7) 


1.06 (4.30) 


1.05 (17.8) 


Fx, hot (10~ n ergs cm -2 s _1 ) 


0.82 (1.79) 


0.89 (6.98) 


0.81 (1.79) 


0.89 (6.81) 



Note. - - Results from simultaneous fits to the EPIC spectra of WR 48a. Tabulated 
quantities are the neutral hydrogen absorption column density (N#), plasma temperature 
(kT), emission measure (EM = j n e TiHedV; see also Fig. [3]) for a reference distance of 
d = 1 kpc, shock ionization age (r = n e t), the absorbed X-ray flux (Fx) in the 0.5 - 10 
keV range followed in parentheses by the unabsorbed value (F x ,hot denotes the higher- 
temperature component, kT 2 ; or kT> 2 keV for CIE Plasma an d NEI Shocks). The 



deriv ed abundances are with respect to the typical WC abundances (Ivan der Hucht et al. 
19861 ) . Errors are the la values from the fits. 



